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ABSTRACT 

A  dynamic  inventory  control  model,  formulated  in  terms  of 
discrete  variable  servomechanism  theory,  is  examined.  The 
behavior  of  the  system  is  first  analyzed  by  digital  computer 
simulation.   Both  random  and  deterministic  inputs  are  considered. 
Two  double  echelon  models  are  then  constructed  and  examined 
analytically  utilizing  transform  theory.  Finally,  the  proper- 
ties of  forecasters  of  demand  are  analyzed  for  the  cases  where 
mean  demand  is  constant  and  linear.   Exponential  smoothing 
techniques  are  compared  with  standard  estimation  procedures. 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA  93943-5101 

TABLE  OF  CONTENTS 

Section  Page 

1.  Introduction  7 

2.  Digital  Computer  Simulation  13 

3.  Double  Echelon  Models  18 

4.  Demand  Forecasting  34 

5.  Conclusions  and  Recommendations  38 

6.  Acknowledgements  40 
Bibliography  41 
Appendix  I -A  42 
Appendix  I-B  43 
Appendix  I-C  44 
Appendix  I-D  47 
Appendix  I-E  51 
Appendix  I-F  55 
Appendix  II-A  59 
Appendix  II-B  60 
Appendix  II-C  61 
Appendix  II-D  62 
Appendix  II-E  63 
Appendix  II-F  64 
Appendix  III -A  65 
Appendix  III-B  68 
Appendix  III-C  72 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

I*  Flow  Graph  Diagram  of  the  Basic  Model  9 

2.  Flow  Graph  Diagram  of  Retail/Wholesale  20 

Double  Echelon  Model 

3.  Block  Diagram  of  Stock  Point/ Inventory 

Control  Point  Double  Echelon  Model  30 

4.  Flow  Graph  Diagram  of  the  Inventory 

Control  Point  31 


1.   INTRODUCTION 

This  paper  will  take  for  its  starting  point  a  dynamic 
inventory  control  model  formulated  in  servomechanism  theory 
by  Vassian  [6]  and  extended  by  Reilly  [5]  and  Zehna  [7].  The 
techniques  of  linear  system  analysis  as  developed  by  Howard  in 
[2]  and  the  exponential  smoothing/forecasting  techniques  of 
Brown  [ 1]  are  heavily  relied  upon.  The  model  as  it  is  described 
in  [7]  will  henceforth  be  referred  to  as  the  "basic  model." 
Briefly,  in  the  basic  model  we  have  demand  as  the  input  and 
inventory  level,  measured  relative  to  some  safety  level,  as  the 
output.  Inventory  is  replenished  using  a  reorder  rule  which  is 
based  on  a  linear  combination  of  past  demand  and  inventory  level 
feedback.  Units  of  inventory  are  assumed  to  be  issued  instan- 
taneously and  a  fixed  lead  time  is  assumed  for  receipt  of  reple- 
nishment orders.  In  discrete  variable  servo  theory  notation  the 
two  equations  which  describe  the  basic  model  are 


(1-1)    I-I.+9.-X     and 
t    t-1    t-l-T    t 


't  mL  Vt-j  +  L     j  t-j 


(1-2)   et 

j-o        j=o 


where    I  ■  on  hand  inventory  relative  to  some  safety  level 
at  time  t. 

6  -  reorder  quantity  computed  at  time  t. 
t 

(A  }  -  a  set  of  constants  which  weight  past  demands. 


{b  }  s  a  set  of  constants  which  weight  past  inventory 

J 

levels* 

Let    D(z)  s   di.  I   >  *>e  c^)e  ^-transform  of  the  sequence  C^vJ 
k«o 

Then,  the  above  equations  become,  in  transform  notation, 

|       |  T+l 
(1-3)  I(z)  -  *I(z)  i  z        9(z)  -  x(z)   and 


(1-4)  9(a)  -A(z)X(z)  +  B(«)I(z>. 

It  is  easily  shown  that 


«T+1 

(1.5)  I(2) ,  m  *<«>  -i  i 

1-1  -  |    B(z) 


Let 

u*    A(z)  -1  ] 
(1-6)   S(z)  «-^ 


X(z) 


I  iT+1 
1-4;  «    B(z) 


be  called  the  transfer  function.  Then. (1-5)  may  be  written 
I(z)  -  S(z)X(z). 

Similarly, 

(1-7)  9(z)  -  CA(z)  +B(z)  3(z)]  X(z) 

«  s'(z)  X(z) 
where   s'(z)  ■  Ca(z)  +  B(z)  S(z)] 

the  basic  model  is  summarized  in  flow  graph  form  in  Figure  1* 
As  Vassian  suggests  [6],  we  will  let 

(i-8)   b(z)  -  iiI"*L- 

,  T+l 
1-z 


:'.: 

i 

I 

.. 


It  follows  then  that  S(z)  will  have  no  finite  poles.  In  that 
case,  the  denominator  of  S(z)  will  be  finite  for  all  finite  values 

of  a. 
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The  following  additional  relationships,  developed  in  [7], 
are  listed  for  ready  reference: 

(1-9)     E(z)   ■  e[i(z)]   -  S(z)  m(z) 

where 

m(z)   -  E[x(z)] 

t+1  a«*T+S 

(1-10)     A*(z)=z  A(z)     iii_->   X(z) 

1-z 

T 

c 

(1-11)  A*  =1  +  )     X     , 

t+l-FT  t+l+T       L       t+l-HT-j     .     We  require  that 

j=o 

t+l-KT 

(1-12)   Efa*  3  -  /     m(j)-*0     as  t-^« 

t-l-l-HT  " 

j-t+1 

When  expected  demand  does  not  vary  with  time  (eCx  ]  =  a)  we  have 

the  following  additional  relationships. 

t 
a  V7     k 

(1-13)  a  ■  a  )      £  X  ,   using  exponential  smoothing, 
t     tmi  t-k 

k«o 


a  x(z) 
(1-14)  a(z)  =   " ■-■V-J-  " 

1  -  0  z 


(1-15)  A*      ■»  (T-:-l)  a  ,  a  forecast  of  demand  during  lead  time 
t-fl-l  |        t 

plus  one  unit  of  time. 
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(1-16)  A*(z> 


T+l 
a  (T+l)    2         X(z) 


1-Pz 


■/ 


(1-17)  A(z) 

a  (T+l)    (l-z) 

T  1-1 
(l-z      )(1-  ?() 

(1-18)   S(z) 

T+l 
a  (T+l)    * 

T+l 
(1   -*1      ) 

* 

1-  pz: 

■ 

l-z 

(1-19)   I(z)   = 


or  <I+1)  ^+1  X(z)    _(1.  J+1)X(Z) 


1-    p'  z. 


1  - 


(1-20)    9(z)   =  s'(z)  X(z)   - 


ot  (T+l)    (l-z),    , 

1  -  i 


X(z) 


Inventing  the  transform  in  (1-20),  we  obtain 


t  t-1 

(1-21)    9     -  a  (T+L)  yV  x         -  a  (T+L)  }'    3     X     ,     +  X 
t  *->  t-k  u  t-l-k       t 


k*o 


k-o 


and 


(1-22)   6     -  a  (T+l)   [x„  -  a     J  +  X     . 
t  t         t-1  t 


Inverting  the  transform  in  (1-19)  yields 


t-T-1  T 

(1-23)  I  -  or  (T+l)  r  P  X        -  Y   Xt_k 

t                      ,  t-T-l-k   , 

k°o  k=o 
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Computing  the  expected  value  of  I  we  obtain 


t-T-1  T 

(1-24)  E  -  a  (T+l)  r   &k  m(t)  -  f  m(c> 

t                Li  Lj 

k»o  k»o 


t-T-1 

r*   k 
«  a?  (T+l)  a\   p  -  (T+l)  a 

k=o 


r  st_T    n 

-  a   (T+l)  a  — — —  "  (T+l)  a   or, 

P  -  1 

t-T 
(1-25)  E  -  -a  (T+l)  3 


Thus,  we  have  described  the  basic  model  and  displayed  the 
resulting  mathematical  relationships.   In  addition,  for  the 
case  of  constant  mean  demand  we  have  derived  explicit  formulae 
for  the  expected  on  hand  inventory  and  reorder  quantities  at 
time  t  when  exponential  smoothing  is  used  to  estimate  the  mean. 
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2.   DIGITAL  COMPUTER  SIMULATION 

The  simulation  presented  here  was  programmed  in  FORTRAN 
60  and  run  on  the  CDC  1604  serial  number  1.  The  programs  used 
are  included  in  Appendix  II  in  the  order  of  discussion  below. 

First,  a  deterministic  simulation  of  the  basic  model  with 
constant  mean  was  made  based  on  the  relation 

t-T-1        T 

E  -  <*(T+1)  Y      Bk  m(t)  -  ^  rcCt) 
k=o  k-o 

letting  at  -  .2,  B  -  .8,  and  T  »  4.  in  the  first  run  m(t)-16 
and  the  results  of  995  iterations  are  shown  in  Appendix  I-A. 
It  is  seen,  that  to  six  decimal  place  accuracy,  the  expected 
inventory  level  reaches  zero  after  83  time  periods.  This  is 
merely  a  simple  verification  of  the  analytical  results  shown 
in  [7]. 

In  the  second  deterministic  run,  we  let  m(t)»t  in  the  basic 
relation.  Since  Reilly  [5]  has  shown  that  first  degree  expo- 
nential smoothing  yields  biased  results  when  the  demand  is 
linear  we  would  expect  this  bias  to  be  reflected  in  the  expected 

inventory.   In  the  case  at  hand  we  have 

t-T-1  T 

Et  -  *(T+1)  £   Bk  m(t-T-l-k)  -  £m(t"k) 
k«o  k=° 

Letting  a   -  .2,  B  ■  .8,  m(t)  -  t  and  T  ■  4, 

t-5  4 

Et  »^(.8)k  (t-5-k)  -]T<t-k) 
k=o  k=o 
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t-5 


t-5 


t-5 


tY'(.8)k-5r(.»k-rk(-»-rt+rk 


k-o 


k=o 


k=o   k=o 


t-4 


.2 


l-<-8) 


.2 


t-4 


1 


+  .8 


t-4  t-5 

(1-.8   )  -  .2(t-4)  .8 

.04 


-5t  -!-  10 


t-4        t-4        t-5        t-5 
-  -5t(.8)    +  45(.8)    -  4t(.8)    +  16(.8)    -35. 


Hence,  lim  E  ■  -35 

t—  t 


We  can  see  from  the  results  shown  in  Appendix  I-B  that  in 
fact  E  ■  -35  after  78  iterations  again  to  six  decimal  place 
accuracy.  So  we  have  verified  the  analytical  conclusions  in 
[5]  and  [7],  i.e.,  that  the  model  is  asymptotically  unbiased 
(no  steady  state  error  in  the  expected  inventory  level)  when- 
ever asymptotically  unbiased  forecasters  are  used  and  asympto- 
tically biased  when  correspondingly  biased  forecasters  are  used. 
However,  as  a  periodic  review  inventory  model,  the  basic  model 
may  require  an  impractical  number  of  periods  to  converge  to  an 
acceptable  level.  Conceding  that  convergence  can  be  ace 
by  adjusting  ex   and  T,  we  must  face  the  fact  that  asymptotic 
unbiasedness  may  not  be  a  practical  measure  of  effectiveness  for 
an  inventory  model. 
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The  second  simulation  to  be  considered  is  stochastic  in 

nature.  Assume  that  demand  is  Poisson,  i.e., 

-Xt  k 

e    Xt 

P[X  -k]  =  - where  X  ■  X  . 

t  k!  t 

For  ease  of  simulation  a  normal  random  number  generator  was 
emp|oyed  taking  advantage  of  the  fact  that  for  large  A.  [3] 


e       \ 

r  k-\  +  .5 

fT           e"^y2     dy 

k! 

inr 

'     k-\  -   .5 

The  1000  random  numbers  for  the  first  two  runs  are  listed  in 
Appendix  III -A,  while  those  used  in  runs  three  and  four  are 
shown  in  Appendices  III-B  and  III-C  respectively. 

In  the  first  run  X  -  16  in  the  basic  relation 

t 
t-T-1  T 

k 


x-«cm>£  p  xt-T-i-k-rxt-k- 


k=o  k«o 

I  thru  I     are  listed  from  right  to  left  in  Appendix  I-C. 
6       1000 

It  is  seen  that  there  is  some  tendency  toward  stability  through 
I-  but  after  that  time  period  we  observe  the  effects  of  the 
random  demand  on  the  inventory  level.   It  is  now  clear  that 
carrying  expected  values  through  the  model  does  not  adequately 
describe  the  true  situation. 

In  the  second  stochastic, run,  all  parameters  and  variables 
were  unchanged  from  the  first  run  except  that  T  was  increased 
to  T  ■  10.  We  see  the  resulting  inventory  values  in  Appendix  I-D. 
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As  expected,  some  increase  in  the  variation  of  I  is  observed, 
i.e.,  the  variance  of  the  random  variable  I  is  a  function  of 
the  lead  time  of  the  basic  modal.   However,  we  also  observe 
that  I  is  relatively  insensitive  to  T  in  the  sense  that  the 
order  of  magnitude  of  the  change  in  T  is  not  carried  over  to 
I  to  an  appreciable  degree. 

In  the  third  stochastic  run,  all  parameters  and  variables 
were  the  same  as  in  the  first  run  except  that  X  =  100.  Natur- 
ally, to  preserve  our  Normal  approximation^  the  Poisson  distri- 
bution variance  was  also  increased  to  100.  Note  that  the  random 
numbers  generated  and  listed  in  Appendix  III-B reflect  this 
change.  Now  we  observe  the  resulting  values  of  I  shown  in 
Appendix  I-E.  Analogously  we  see  that  the  variance  of  I 
depends  on  the  variance  of  X  but  is  relatively  insensitive. 
The  order  of  magnitude  of  the  increase  in  X  is  not  reflected 
in  I  . 

In  the  final  run  X  was  allowed  to  vary  with  time,  speci- 
fically X  -  16  +  t.  From  the  analytical  development  we  are 
led  to  expect  bias  since  our  forecast  is  biased  for  non-con- 
stant mean  demand.   In  Appendix  III-C  observe  the  new  set  of 
random  numbers  generated  with  mean  X  =  16  :  t  and  variance 
X  «  16  +  t.  The  expected  bias  is  apparent  from  an  inspection 
of  the  Appendix  I-F.  Most  values  of  It  are  negative.  Again 
observe  the  relative  insensitivity  of  I  to  increasing  the 

variance  of  the  demand. 

In  summary,  the  computer  simulations  have  confirmed  the 
theoretical  results  as  expected.,/  However,  in  doing  so,  some 
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limitations  of  the  basic  model  have  been  suggested.  Asymptotic 
unbiasedness  as  a  criterion  appears  to  be  weak  and  probably 
impractical  even  when  expected  values  of  the  random  variables 
X  and  It  are  considered.  When  {x  }  is  a  sequence  of  realizations 
of  a  random  variable  from  a  probability  distribution,  the  sequence 
{i  }  certainly  does  not  appear  to  approach  zero,  Vassian's  [6] 
apparent  claim  to  the  contrary.   Zehna  [7]  clarified  the  theory 
in  this  area  but  the  results  of  the  simulation  dramatically 
emphasize  the  inadequacy  of  the  basic  model  when  only  expectations 
are  considered.  The  choice  of  a  safety  level,  or  other  managerial 
decisions,  cannot  be  made  on  the  basis  of  E[l  ]  alone  but  Var[l  ] 
must  be  accounted  for  also. 
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3.  DOUBLE  ECHELON  MODELS 

It  is  possible  tie  formulate  multiple  echelon  models  in 
servo  theory  along  the  lines  suggested  by  Howard  [2]*   It  is 
desired  to  test  the  basic  model  of  the  previous  sections  in 
echelon.  First,  we  construct  a  double  echelon  model  using 
three  basic  models,  two  at  the  retail  level  and  one  at  the 
wholesale  supporting  stock  point  position.  Assume  that  con- 
sumer demand  with  constant  mean  values  a  and  a  is  the  system 
input  at  the  retail  level.  In  addition  to  the  inventory  level 
outputs  at  the  retail  stock  points  we  have  the  reorder  quantities 
9  (t)  and  9  (t)  serving  as  demand  inputs  to  the  wholesale  stock 
point.  There  is  a  lag  time  which  is  constant  for  each  retail 
stock  point  of  T.  and  T  respectively  between  the  time  of  issue 
by  the  wholesale  stock  point  and  receipt  by  the  retail  points. 
Assume  that  the  wholesale  stock  point  has  a  sufficiently  high 
safety  level  so  that  it  is  always  able  to  fill  orders  from  the 
retail  activities.  There  are  to  be  no  interactions  between  the 
retail  stock  points*  All  three  stock  points  are  to  use  a  single 
exponential  smoothing  forecaster  which  is  of  course  asymptotically 

unbiased  for  constant  mean  demand.  The  lead  time,  T  .  for  the 

3 

wholesale  stock  point  is  the  delivery  time  from  the  manufacturer  . 
Assume  that  all  three  stock  points  use  the  same  value  of  the 
smoothing  constant  ot.     Hence  this  double  echelon  model  describes; 
in  a  simple  way,  a  common  situation  involving  two  eel  elorts  of 
inventory.  The  use  of  a  common  inventory  control  doctrine  and 
the  assumption  of  no  interactions  at  the  retail  level  are  considered 
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to  be  reasonable  In  the  light  of  the  author's  experience*  It 
is  then  of  interest  to  determine  the  possible  effects  of  such 
an  overall  doctrine  in  terms  of  this  model. 

In  the  framework  of  the  multiple  input-output  system  in 
[4]  we  have  a  system  with  two  inputs  X  (z)  ,  X  (z)  and  three 
outputs  I,(z),  I~(z)  and  I  (z) .   In  matrix  notation  this  becomes 


[i^z)  I  (z)  I  <z)]  =  [xi(z)yX2(z)] 


T 
h  (z) 
11 


T 

h  (z) 
22 


T 
h     (z) 

13 

T 

h     (z) 
23 

where,  for  the  present  case, 

Tl  +1 

z  cv(Ti    :  i)        1-z 


IX  :■  1 


h     (z)   -  Sfz)   = 

11  1 


1  -  0  z 


1-z 


T2  +  1 


•£  Z 

and     h     (z)   ■  S   (z)   ■  — 

22  2  1  -  0  z 


"<T2+1>_   1-zV1 


1-z 


Hence,   I  (z)   -  S   (z)  X  (z)  and 


I  (z)   -  So(z)   X9(z)      . 
2  2  2 


This  double  echelon  model  is  graphed  in  Figure  2.     Now 

I   (z)    -  tJ„(z)    [6   (z)  +  9  (z)] 
3  ->  1  2 


since  the  input  at  the  wholesale  stocking  point  is 


e^z)  +  e  (z) 
1     2 


19 


Us 

2n 


0= 


J 


H 


F  1 

It  "^ 

ft 

t 

j 

., 

cvl 

Or 

CD 


3 

CO 


►4 
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From  (1-7)  we  have 


I  (z)   -  S   (z)   {[A(z)  +  B.(z)   Sfz)]  X,(z)  +  U  <z)+B  (z)S   (z)] 
3  3  1  111  222 


X2(z>} 


Hence > 


h     (z)    »  S.Cz)   UAz)  +  B,(z)   S.(z)]  and 

13  3  1  1  A 


n     (z)   -  S   (z)   [A  (z)   +  B   (z)   S0(z)]    . 
23  3  2  2  2 


Recall  that  we  arc  assumimj 


m,(t)    =  E[x,(t)]   =  a,   and  m  (t)   -  e[x  (t)]  »  a     for  t  *  0. 
1  112  2  2 


Hence , 


i  2 

m,(z)    = and  m  (z)   ■  • 


1-z 


1-z 


Consider  now 


t-1 


9.(t)   -  «  (Ti+1)  I  3     X^k  -  or  (T.+l)  ^  P     Xt.l-k+ 


X(t) 


k=o 


so  that 


e^z) 


orCij+1)    (1-z) 

i +   1 

1  -  3  z 


k=o 


X.(z) 

i 


t-1 


for  i  -  1,  2. 

t 

E[9   (t)]  -  a(T  +1)   a.  )   &     .  a(T  +1)   .    VP     +  a 


k«o 


k=o 
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t  t-1 

(9  -l)  (9   -l) 


E[e  (t)]  =  or(T  +1)  a^  -  or(T  +1)  a^  +  a, 


9-1        i     4   9-1 


t-1 


all-CT^l)  Ofc-l)  +  (T±+l)  O  "  -1)] 


a  Cl-(T  +1)  9  +  (T  +  1)  9    ] 
i     i  i 


a  [1+(T  +1)  (9*"  -  9*)] 
i     i 


t-1 

»  a  [1+(T  +1)  9    (1  -  9)] 
i     i 


t-1 

a  [1+(T  +1)  <*  9    ]  . 


This  demonstrates  that  the  expected  demand  at  the  wholesale 
stock  point,  which  is  the  sun  of  the  order  quantities  from  the 
retail  stock  points  [9  (t)  +  9  (t)  ],  is  not  a  constant  despite 

the  constant  mean  demand  inputs,  m  (t)  and  m  (t)  ,  to  the  system 

1        2 

as  a  whole.  It  is  clear  then  that  the  inventory  control  doctrine 
at  one  echelon  does  affect  the  stability  of  the  system.  In  this 
case,  a  constant  mean  demand  is  transmitted  as  an  exponential 
function  because  of  the  exponential  smoothing  forecasting  system 
employed. 

Now  it  is  of  interest  to  examine  the  behavior  of  the  whole- 
sale inventory  level  under  these  conditions.  Consider 
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S   (z)    [A(z)   +  B_(z)   S   («)]   X(z) 
3  1  111 


which  is  the  first  component  of  I  (z)   shown  above.     By  assumption 


s3(^)  - 


<^(T3-H)    z 


T3+l 


1-z 


T3+l 


1  -  3  z 


1-z 


[A^z)  +  Bx<z)   S^z)]  X1(z) 


T,+l 


aCT^l)    (1-z)  z  l    orCT^lXl-z) 


T  +1 
(1-z  l    )(l-3z) 


T,+l 
(1-z  1     )    (1-3  z) 


+  1 


1 
1-z 


o/CTj+1)   (1-z)  +  (1-3  z) 
1  -  3  z 


a 
.        1, 

1-z 


Hence, 


S   (z)[A1(z)   +  B,(z)   S.(z)]  X  (z)   =  S   (z)   s'(z)   X  (z) 
■al  x  1  1  311 


T3+l 
a(T3+D    z 

T3+r 

1-z 

^^(l-zJ+tl-Pz) 

1  -  3  z 

— 

1-z 

1  -  3  z 

-                                            * — 

a 
1 

1-z 
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Now, 


lim 
z-^  1 


T  +1 
^(T3-H)z  3 

1  -  3  z 


1-z 


T3+l 


1-z 


-  0  from  [7,  p.  15], 


and  obviously 

<*(Tril)  (1-z)  +  (1-3  z) 


lim 

z  ->  1 


1  -  3  z 


1  . 


Therefore, 


lim  (1-z)  S  (z)  -  0 
z   1        3 


and  since 


E  (z)=So(z)  Cs'(z)  X  (z) 
3     3      11 


S^(z)  X2(z)]  , 


it  follows  that 

lim   (1-z)  E  (z)  -  0 

z  -:   1        3 

Based  on  a  well  known  theorem  |4,  p.  218],  the  Result  above 

implies  that  E  — ^0  as  t->°°. 
t 

At  first  this  may  seem  to  be  unexpected  until  one  realizes 
that  asymptotic  unbiasedness  is  not  a  necessary  condition  for 
the  steady  state  error  to  vanish.  Recall  that,  in  transfer 
notation, 

(l-z)E  (z)   «  S  (z)s'(z)  a,  +  S  (z)s'(z)  a  . 

3       3    l  1    3    ?     l 


Thus  with  constant  mean  demand,  the  most  that  is  required  for 

asymptotic  unbiasedness  is  that  lim   S  (z)  *=  0  and  that 

z  ->1   3 

lim  S j  (s)  be  a  constant.   Since  S  (z)  depends  only  on  A(z) 

z   1  3 
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and  B(z) ,  it  is  relatively  easy  to  choose  A(z)  and  B(z)  to  yield 

the  limit  zero.  S.(z)  here  is  a  linear  combination  of  the 

asymptotic  unbiased  retail  level  transfer  function  and  hence  will 

be  a  constant  in  the  limit  under  these  conditions*  Therefore. 

the  results  above  are  indeed  natural  under  the  circumstances*  If 

initial  expected  demand  is  not  a  constant  (m(z)  ^  a  /(L-z)  then 

S  (z)->0  and  S  (z)->c  for  some  constant  c  may  not  be  sufficient 
3  i 

conditions  for  no  steady  state  error*  Again  this  emphasizes  the 
weakness  of  asymptotic  unbiasedness  as  a  criterion* 

Inventory  managers  have  a  keen  interest  in  the  number  of 
units  in  the  pipeline  (units  issued  from  one  stock  point  but 
not  reflected  on  the  records  of  another).  Let  us  see  what  this 
double  echelon  model  reveals  about  the  problem*  Remember  it  was 
assumed  that  issues  are  made  instantaneously  by  the  wholesale 
stock  point  but  T  +1  and  T  +1  time  periods  elapse  before  the 
inventory  is  recorded  at  the  respective  retail  points.  We  are 
interested  in  the  expected  number  of  units  in  the  double  echelon 
system  which  are  not  shown  on  the  records  of  any  stock  point  as 
inventory. 

The  total  number  of  units  on  order  at  time  t  by  a  retail 
stock  point  can  be  written 

T       T         t-j 

k 
8  X 

-k 


j=l      j=l        k»o 


Hence , 

T    t-j 

B[Q  ]  -  <*  (T  +  1)  V       £  Bk  a 
j-1    k-o 
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T  t-j 

B[  Ql]  =  a(T-l-l)    a  X1  V    ak 

L  L    p 

j=l  k=o 


a-(T+l)    a 


T  t-j+1 


La 


1-1  P  -  1 


T     t-j+1  ^    t-j+1 


=  -a(T+l)  V(P  -1)    =  aT(T+l)    -a(T+l)  /P 

j-l 


j-l 


t+1                                    t-T+2 
=  aT(T+l)-a(T+l)   — —     +  a(T+l)  — — 

p   -   1  P   -   1 


aT(T+l)    +  a 


(T+1) 


a 


t+1     t-T+2 
P        -P 


t-T+2 
a(T+l)    P  T-l 

aT(T+l)   + (P       -1) 


a 


=  aT(T+l)    - 


t-T+2         T-l 
a  (T+1)    P  (1-p       ) 


<x 


Then,    lira     e[q  ]   -  aT(T+l) 


for  any  retail  point  with  e[x  ]  ■  a 


and   so 


lim     E[Q  ]   =  a  T   (T     +   1)   +  a  T   (T     +  1) 
t  111  222 


t      .-co 


for  the  present  double  echelon  model. 
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This  is  a  natural  result  and  the  steady  state  condition  is 
a  direct  consequence  once  again  of  the  asymptotically  unbiased 
forecasting.  Note  that  the  expected  amount  in  the  pipeline  is 
a  function  of  the  square  of.  l;:o  1  time  w  ch  re -emphasizes  the 
possible  savings  resulting  from  a  reduction  in  the  lead  time  in 
any  real  situation.  Then  it  appears  that  models  of  this  type 
are  useful  for  further  study  of  pipeline  inventory. 

An  inventory  control  point  (ICP)  is  an  office  whose  function 
is  to  manage  a  so-called  commodity  area  which  usually  consists  of 
a  large  number  of  stock  items.  A  substantial  portion  of  its 
effort  involves  the  maintenance  of  optimal  inventory  levels  at 
a  system  of  stock  points.  The  stock  points  report  changes  in 
their  inventory  levels  and  forward  replenishment  orders  to  the 
ICP  periodically.  The  ICP  takes  centralized  procurement  action. 
One  point  of  great  interest  is  the  correlation  between  inventory 
at  the  ICP  and  the  inventory  as  it  really  exists  at  the  stock 
points.   It  is  possible  for  incompatible  inventory  control  doc- 
trines at  the  various  stock  points  as  well  as  different  lead 
times,  to  generate  a  record  of  demand  and  inventory  at  the  ICP 
which  differs  greatly  from  the  actual  demand  input  to  the  system 
and  the  inventory  at  the  stock  >oints,  A  very  simple  double 
echelon  model  which  operates  on  the  basis  given  above  can  be 
analyzed. 

Assume  there  are  two  stock  points  operated  as  basic  models 
with  constant  mean  demand.  Again  the  mean  demands  a.  and  a  are 

not  necessarily  the  same  nor  are  the  lead  times  T  and  T  •  There 
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are  no  interactions  between  the  stock  points.  The  reorder 

quantities  generated  at  the  two  stock  points,  0  (t)  and  G  (t) , 

1        2 

are  transmitted  to  an  inventory  control  point,  whose  operation 

differs  from  that  of  the  basic  model.   It  is  assumed  that, 

unlike  the  stock  points,  the  control  point  has  no  physical 

inventory.  It  merely  keeps  a  record  of  system  inventory  carried 

at  the  stock  points  based  on  the  inputs  it  receives.  Besides 

the  two  order  quantities,  the  ICP  receive;;  a  report  of  the 

inventory  level  at  each  stock  point  at  each  time  period.  Material 

ordered  by  the  ICP  is  delivered  direct  to  the  stock  points  so 

that  the  lead  time  T  is  a  part  of  T  and  T  .  The  control  point 

does  not  feed  back  its  recorded  inventory  level  but  utilizes  the 

values  of  I„(t)  and  I  (t)  to  determine  order  quantities  6  (t) 
12  3 

and  update  I_(t) .  The  basic  equations  which  describe  the  oper- 
ation  of  the  ICP  are  therefore 


I,(t)  =  In(t-1)  +  I  (t-1)  :  8  (t-l-T  )  -  [9  (t)  +  9  (t)] 
3      1        2        3      3      12 


t  % 

e3(t)  «£a.  [e^t-j)  +  e2(t-j)]  +^b.  [i^t-j)  +  i2(t-j)] 

j-o  j=o 


and   in  transform  notation 
T3+l 

i  (z)  =  z        e  (2)  -  e  (z)  -  e  (z>  +  z  [i  (z)  +  i  (z)] 

3  3  12  1  2 


9   (z)   =  A   (z)    [9   (z)   +  9   (z)]  +  B   (z)    [I   (z)   +1   (z)] 
3  3  12  3  1  2 
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The  flow  graph  for  this  double  echelon  model  is  shown  in 

Figure  3.  The  inventory  control      ';.  system  is  graphed  in 

Figure  4. 

Now  let  us  examine  the  behavior  of  the  ICF  inventory  record 

in  relation  to  the  actual  inventories* 

T3+l     f 

I,(z)    =  z  [l,(z>+I  (z)]  +  z  J  A  (z)[9i(z)+e  (z)] 

3  1  2  )    3  1  2 


+  b(z)[i  <«>+i  (z)]y-  [e,<»)+e  (z)] 

12  2 

T3+l 

z[s   (z)X,(z)   +  S   (z)X  (z)]  +z        A  (z)[s'(z)X  (z)+s'(z)] 
ll  22  3112 


+  B(z)[s   (z)X  (z)   +  S0(z)X0(z)](-[s'(z)X  (z)+s'(z)X  (z)] 
11  22Jll  22 

T3+l 

[zS   iz)+z         A.(z)s'(z)   +B   (z)S,(z)-s'(z)]  X  (z) 
1  3         i  3         1  1  1 

Y1 

+  CzS   (z>  +  z  *      A  (z>S/(z)   +  B(z)S   (z)-s'(z)]  X  (z) 
2  3         2  2  £  2 


S     (z>X  (z)  +  S     (z)X  (z) 
13         1  23         2 


where 

T  +1 

S     (z)    «=  zS   (z)   +  z         A    ^)s'(z)   +  B(z)S   (z)    -  s'(z)  and 

13  1  3  1  11 

T  +1 

S     (z)   »  zS   (z)  +  z         A  (z)s'(z)   +  B(z)S   (z)    -  s'(z) 
23  2  3         2  2  2 
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Since 


lim   S  (z)  -  0,   lim  A  (z)  -  1, 
z  -*1   1         z  ^1  3 


11m  s'(z)  -  1,   lim   B(z)  = — 

z  >1  1  z  >1         T^l    > 

it  follows  that 

lim   S  (z)  =  1  -  1  -  0. 
z  -^1   13 

Similarly, 

lim    S  (z)  -  0. 
z  -*1     23 


But, 


E  (z)  -  S  (z)  m  (z)  +  S   (z)  m  (z) 
3      13     1      23     2 


al  a 

-  S  (z)  +  S   (z) 


13    1-z     23    i-2 

and 

(1-z)  E  (z)  -  a,  Slo(z)  +  a  S  (z). 
3      1  13      2  23 

Therefore, 

lim   (1-z)  E  (z)  -  0, 
z-^1         3 

which  implies  that 

lim   E  (t)  ■  0. 

t  — >   CO      3 

Again  there  is  no  steady  state  error.  Basically  the  results  here 
parallel  those  in  the  first  double  echelon  model.  In  effect, 
major  difficulties  can  exist  in  tlse  type  of  system  iscussed 
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above  but  the  "no  steady  state  error"  criterion  does  not  really 
measure  the  behavior  of  the  output  function  adequately.  It 
appears  feasible,  however,  to  formulate  multiple  echelon  servo- 
mechanism  inventory  models  which  include  ICP  type  operations  in 
order  to  gain  further  insight  into  the  operation  of  such  a 
complex  ;     u 
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4.   DEMAND  FORECASTING 

Reilly  [5]  and  Zehna  [7]  apply  exponential  smoothing  to 
define  forecastors  in  the  cases  where  mean  demand  is  constant 
and  also  when  it  is  linear.  Moreover,  in  terms  of  the  basic 
model  presented  here,  Zehna  [7]  has  demonstrated  that  the 
expected  inventory  level  converges  to  zero  whenever  any  asymp- 
totically unbiased  forecastor  is  used.  For  example,  in  the 
constant  and  linear  models,  if  least  squares  estimators  are 
used  to  estimate  the  parameters  and  these  are  substituted  in 
the  forecastor,  then,  being  unbiased,  the  forecastor  is,  a 
fortiori,  asymptotically  unbiased  and  the  inventory  level  will 
converge  to  zero.  The  justification,  then  for  using  exponential 
smoothing  lies  in  the  advantages  which  relate  to  digital  com- 
puter storage  and  mathematical  tractability  in  the  servomechanism 
model.   See  Brown  [ l]  for  a  complete  treatment  of  the  advantages 
of  exponential  smoothing.  However,  these  advantages  are  external 
to  the  present  requirement  of  unbiasedness  in  the  limit. 

Under  the  model  m(t)  ■  bt,  Zehna  [7],  following  Reilly  [5], 

uses  b  =  -3-  [s^(X>  -  S   (X)]  where  S   (X)  represents  the 
t    P    t       t  t 

A 

sequence  X  doubly  exponentially  smoothed.   It  follows  that  b 
t  t 

is  an  asymptotically  unbiased  estimator  for  b  so  1     tjhe  '. 
caster        t+T+1 


A*    -l        V 
t+1«   J-t+1 

is  an  asymptotically  unbiased  estimator  for  total  demand  over 
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the  next  lead  time*  For  this  model,  another  possibility  suggests 
itself  quite  naturally*  Since, 

E[S  (x)]  -  S  <S)  -  bt  -  — 

t       t  a 


E[s(x)  -  x(t)]  =  -  -3- 

t  G, 


and 


a 


11m   E[  —  (X(t)  -  St(x))]  =  b 


Hence,  the  estimator 

C  or 


$  Ex(t)  -  st(x)] 

is  also  asymptotically  unbiased  for  the  parameter  b.  Using  this 


to  define  a  forecaster 
t+T+1 

A*  rT  tj 

t+T+l    ^  tJ 

j=t+l 


it  follows  that 


e[a*       ]  »  e 

t+T+1 


t+T+i  t+T+i        t+r+1 

j»t+l       j=t+l      j=t+l 


so  that  A*      defines  an  asymptotically  unbiased  estimator  for 
t+l+T 

demand*  Again  by  Zehna  [7],   lim  E  ■  0  when  this  forecaster 

t  -*00   t 

is  used.  Since  this  basic  result  is  the  same  as  that  obtained 

when  double  exponential  smoothing  is  used,  other  means  will  have 

to  be  used  for  comparing  the  two. 

First  consider, 

E|^[st(*)-sCt23(x>]]--2  [st(l)-st2  (5)] 
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^  ,tx         bP  r         t   .  t-1. 


t  a  <*L     p  J  „ 


Pa  2 


t+1  t 

a     fbp_  .  b§ t(t+l)b  otp 


p         a 


2  t 

t       t(t+l)b  a  P 

o    b    -    bP       - 


2  P 


On  the  other  hand 


or 


fi        a 


E  J  —  [X(t)    -  S   (x)] .,  =     a-  [bt  -  S  (§)] 


p  t    J    p 


t 


(4-2) 


—  fbt   -   I   (1-pV  ^  [1-pV-t)    -  tP        ] 
P    I 


2L(bt  -  bt  +  btp*+  —  [l-p'u-tf-tp'"  ] 

P    I 


t 


t        *        „t-l. 


f-jbtP    -i-bP[i-3C+  tP     -  tp'"1] 


2L  c  M         t^  +  bP  ^t-lj  +        t] 
p         or  or 


36 


t 
i/-2  [x(t)  -  s(x)]l  -   -f  [  ^  O-P*)  •  *tjL  O-D+  bt3*  ] 


p  tWJ  \  P  L     or 


P        or  P        cr  a 


b  -  bp     . 


Obviously  as  t-lf0  (4-1)   converges  to  b  more  rapidly  than  (4-2) 
does  and  so  is  preferable* 
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5.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  conclusions  reached  here  are  that  expected 
inventory  levels  and  asymptotic  unbiasedness  appear  to  be 
rather  weak  criteria  for  the  models  considered  in  this  paper* 
It  would  be  highly  desirable  to  be  able  to  relate  costs  and 
effectiveness  to  these  models.  The  inventory  variance  should 
be  the  basis  for  any  discussion  of  effectiveness  (expected  per-  - 
centage  of  demands  issued)  and  the  quality  of  the  estimator. 
Average  cost  depends  on  inventory  carrying  costs,  order  costs 
and  shortage  costs  mainly.  It  is  important  that  the  decision 
criteria  be  related  to  these  costs. 

Many  possible  courses  of  action  involving  the  ramifications 
of  the  basic  model  exist.  It  is  recommended  that  the  following 
two  extensions  be  attempted  first  to  determine  whether  or  not 
this  type  of  model  has  sufficient  power  to  be  useful: 

(1)  Calculate  the  Var[l  ]  for  several  types  of  random 
demand  and  forecasting  schemes.  Explore  the  methods  of  deter- 
mining a  safety  level  to  provide  a  given  degree  of  protection 
against  shortages. 

(2)  Attempt  to  define  a  new  criterion  for  acceptable 
inventory  level  behavior  which  can  be  related  to  the  commonly 
considered  costs  discussed  above.  Some  insight  into  the  pro- 
blem might  be  gained  by  relating  the  lack  of  steady  state  error 
criterion  to  costs.  Specifically,  cost  would  seem  to  depend  on 
the  amount  the  inventory  level  varies  from  the  norm,  the  relative 
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values  of  the  cost  components,  the  number  of  times  I  changes 
direction  and  the  rate  of  convergence  of  E  • 

For  this  model  to  be  useful  it  is  necessary  that  the 
reorder  policy  be  optimal  in  the  sense  of  minimum  cost.  After 
a  minimum  cost  criterion  is  applied  to  the  model,  the  stochastic 
process  {i  ,  t  ■  0,1,2, •••}  must  be  intimately  understood. 
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_____    -.000078   -.000063  -.000050  -.000040  -.000032 

"-.0003  26   -.000021  -.000016  -.000013  -.000011 

-.000308   -.000007  -.000005  -.030004  -.000003 

-.000003   -.000002  -.000002  -.000001  -.000001 

_________  -.000001   -.000001  -.000001  -.000000  -.000000 

-.000000  ~-. 000000  -.000000  -.030000  -.000000 

-.000300   -.000000  -.CCOOOO  -.000000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-."OOOOOO-  -.000000  -.000000  -.030000"  "-.000000 

-.000000   -.000000  -.000000  -.030000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-.000000   -.000000-  -.COOOOO  -.000000  -.000000 

~- . 0003  oo -~ - .000000  -. 0000  00  - . 03  0000  " "  - . 000000 

-.000000   -.000000  -.000000  -.030000  -.000000 

-.000000   -.000000  -.000000  -.030000  -.000000 

-.000300   -.000000  -.000000  -.000000  -.000000 

~~ ~      -.0003  00   -.000000  -.0000  00"  -.03  0000  "-.000000 

-.000300   -.000000  -.000000  -.030000  -.000000 

-.000000   -.000000  -.COOOOO  -.000000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

~-\  0003  0 0 ~~  -  . 000000  -.0000  00  "  -  .03  0000  ""-  .  000000 

-.000000   -.000000  -.COOOOO  -.030000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-.000300   -.000000  -.000000  -.030000  -.000000 

"-.000300   -.000000  -.000000  -.ooooon  -.000000 

-.000000   -.000000  -.003000  -.030000  -.000000 

-.000000   -.000000  -.COOOOO  -.000000  -.000000 

-.000000   -.000000  -.000000  -.030000  -.000000 

-.000300   -.000000  -.000000  -.030000  -.000000 

-.000300   -.000000  -.COOOOO  -.000000  -.000000 

-.000000   -.000000  -.COOOOO  -.000000  -.000000 

-.000300   -.000000  -.COOOOO  -.000000  -.030000 

-.0003  00   - . 0 30000  -  . COOOOO  ~-. 03 0000  - . 000000 

-.000000   -.000000  -.COOOOO  -.000000  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-.000000   -.000000  -.000000  -.030000  -.000000 

-.000300   -.000000  -.COOOOO  -.030000  ""-.000000 

-.000000   -.000000  -.000000  -.030C00  -.000000 

-.000000   -.000000  -.003000  -.OOOCOO  -.000000 

-.0003  00   -.000000  -.0000  00  -.03  0000  -.000000 

"-.000000 ""  -.000000""  -.COOOOO  ~~ =7000000  -.000000 

-.000000   -.000000  -.COOOOO  -.OOOCOO-  -.000000 

-.000000   -.000000  -.COOOOO  -.OOOCOO  -.000000 

-.000300   -.000000  -.000000  -.030000  -.CCOOOO 

-.000000   -.000000  -.COOOOO  -.030000  -.000000 

-.000000   -.000000  -.COOOOO  -.OOOCOO  -.000000 

-.000000   -.000000  -.000000  -.000000  -.000000 

-.000300. __-. OOOCOO -.000000 -.-030000 -.000000 

-.000300   -.000000  -.COOOOO  -.030000  -.000000 

-.000000  .  -.000000  -.COOOOO  -.000000  -.000000 

-.0003  00   -.OOOCOO  -.000000  -.OOOCOO  -.OOOCOO 

__-. 000300   -.000000--. 000000  --.OOOCOO-  -.000000- 

-.000000   -.0     iQ  ~. COOOOO  -.OOOCOO  -.000000 

-.000000   -.000000  -.000000  -.OOOCOO  -.000000 

-.000300   -.000000  -.000000  -.000(00  -.000000 

-.0003  00  -.ooooon  -.nnoono  -.ooorno  -.onoono 
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APPENDIX  1-3 


.  MEAN  VALUES 

-19.000000    -22.200000    -24 
-29.757120    -30.805696    -31 

760000    -26.808000    -28.446400 
.644557    -32.315645    -32.852516 

-33.262013    -33.625610    -33 
-34.437050    -34.54964Q    -34 
-34.815533    -34.852426    -34 

900488    -34.120391    -34.296313 
.639712    -34.711770    -34.769416 
.881941    -34.905553    -34.924442 

-34.939554    -34.951643    -34 
-34.980193    -34.984154    -34 
-34.993510    -34.994808    -34 

.961314    -34.969052    -34.975241 
987324    -34.989859    -34.991887 
995846    -34.996677    -34.997342 

-34.997873    -34.998^99    -34 
-34.999303    -34.999442    -34, 
-34.999772    -34.999817    -34 

998639    -34.998911    -34.999129 
999554    -34.999643    -34.999715 
999854    -34.999883"    -34.999906 

-34.999925    -34.99994q    -34, 
-34.999975    -34.999980    -34, 
-34.999992    -34.999994    -34. 

999952    -34.999962    -34.999969 
999984    -34.999987    -34.999990 
999995    -34.999996    -34.999997 

-34.999997    -34.999998    -34. 
-34.999999    .34,999999    _34, 
-35.000000    -35.0000  0  0    -35 

999998  -34.999999    -34.999999 

999999  -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.00000U    -35.000000    -35. 
-35.000000    -35,000000    -35, 
-35.000000    -35.000000    -35. 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35. 
-35.000000    -35,000000    -35. 
-   . 0  0  0  0  C  0    -35.000000    -35. 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35,000000    -35. 
-35.000000    -35.000000   ■ -35. 
-35.000000    -35.000000    -35. 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35 
-35.000000   .-35.000000    -35. 
-35.000000    -35.000000    -35. 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35. 
-35.000000    -35.000000    -35. 
-35.000000    -35.000000    -35. 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35. 
-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 

000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.Q000G0    -35.000000    -35 

000000    -35.000000   "-35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 
000000    -35.000000    -35.000000 

-35.Q0OO0O    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 

000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35,000000    -35 

.000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35 
-35.000000    -35.000000    -35 
-35.000000    -35,000000    -35 

.000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 
.000000    -35.000000    -35.000000 

-35.000000    -35.000000    -35 
-35.000000    -35S.O0OUOO    -35 
-35.000000    -35.000000    -35 

.000000    -35.000000    -35.000000 
OOQ0O0    -35. 000000    -35,000000 
000000   .-35.000000    -35.000000 
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APPENDIX   I-C 
INVENTORY    ON    HAND 


-40. GOOD  00-42. 800000- 32. 840000- 17.87 2 000-29. 897600' 

-40.7180  80-29.5  7  4464-28.4  59571-36. 16  7657-23.934126 

-7.147300    -5.717340      2.G25728    12.220582    11.576466 

4.461173      7.568938      5.655150       5.324  120-1.740704 

-4.192563    11.045950       6.C36760    -5.090592    -9.472474 

-21.177979-12.942383-11.953907       8.236875    18.589500 

18.671500       2.737280    -8.410176-19.728141-15.5825  13 

_-8. 266010    -9.412808       2.469753    15.775303    15.020642 

6.016514     11.013211    -1.989431       2.408455     13.126764 

4.7014 1 1    -6.8  38871-10.4710  97-24.176877-26.341502 

-2.873202      5.301439    12.041151       -.767079       3.186337 

—^-7. 4509 3 1-1 0.360745- 19. 4885 96-- 3. 59 0  877-    1  .927  299 

7.141839      -.286529    17.570777    12.856622      9.085297 

6.868238       3.294590      -.364328       -.691462    -4.553170 

12.557464  .245971    -5.003223    -9.402578    -1.922063 

-■=--12.337650- -3.070120    12.543904-    -.164877  —  5.268099 

15.614479    11.291533    11.833266    11.866613       4.293291 

2.6346  32    -6.492294-20. 1938  35-14.15  5068-13.724055 

-21.179244    -9.943395       9.245284    23.396227       8.116982 

-1.306415  --6. 445132-21. 956105-27. 564884    -6.651907 

5.278474    -1.177221-14.141777       5.886579    -1.690737 

-3.352590       6.717928    15.374343       6.499474       5.799579 

21.239663    10.791731     13.633385    11.906708    12.325366 

3.660293    11.123234       3.5G2587      8.002070     10.001656 

12.601325    -1.518940       1.584348       4.067878    14.254303 

-3. 1965  58    -6.757246-25.0057  97-23.8046  38-53.843  710 

-15.274968       4.180026       6.744020       3.995216       7.796173 

6.436938-5.349551       5.679641  .    6.743712       1.194970 

-7.044024-14.235219  -4.988175-14.390540-10.912432 
1.870054  -2.903957  -5.923165  -7.933532  -8.950826 
-.360661     17.911471       -.270323       2.183342       8.346673 

=8.322661-14.658129    -8.326503-10.86  1203    -1.488962 

3.408330-1.472936  3.021651  3.417321-12.666143 
-7.532914  4.973668  2.178935  13.143148  20.114518 
17.291515  .633292    -9.893367   -9.714693    -6.571755 

__-3. 057404    -9.245923     .1.203262         .962609    11.970087 
23.576070    19.66j356       3.928685   -2.857052-12.285642 
-21.4285  13-11.742811       2.805751    24.644601     16.715681 
22.772545    20.018036       3.414429   -6.468457-10.774766 

_=_17.419813_-7. 935850    -3.348680 .92  1  056.-JO.  336845. 

-5.530524    -6.024419    -8.019535       2.384372       4.907497 

13.725998       8.980798      -.615361    -8.692289    -1.353831 

-11.083065-16.266452      -.413162      9.269471       8.215577 

15.772461_.l .617969-14. 105625-24. 684500-J 1 . 947600 

-11.558080-13.246464  .802329       8.042263       4.233810 

1.987048    13.989639-1.608289       6.313369  .450695 

-1.239444    -8.991555    -4.593244       3.325405    18.260324 

_23. 208259._11.766607      8.213286    -6.229371-10.983497. 

-8.986798      -.389438       9.088449    10.870760       7.096608 

4.077286   3.861829  -5.710537-15.768430-20.814744 

-7.451795    -6.3614,36       2.910851    -4.67  1319-11.937055 

_=J 4. 949644       -.359715"  -4.087772_    4  .  52  9782___1 1  .  223826 

6.179061       2.143248    -4.285401    -5.828321    -6.462657 

-9.370125-10.496100    10.603120    17.032496    31.265997 

26.012797    19.010233       8.403190      6.926552-12.258758 

_=22. 0073  07-12. 20 560 5- 19. 76448 4-24. 61  15  87-15.289270 

-12.031416-13.225133       3.419894         .335915       5.268732 

-2.385314       1.091989       2.273591     10.618873-8.704902 

2.836078       4.868863       3.C95090       -.123928       2.900858 

9.520686  ....  ..616549    -1.506761.,   6.994591 8.195673. 

8.156538    11.325231       7.060185-6.35  1852         .318518 

1.454315       6.763852    -5.388919   -3.911135    -8.123908 

-12.903126-22.122501     -7.298001       4.161599       1.729279 

=5.416576 8.866739  _     .693391 ..-5.  64  5287_-9. 316230 

-6.052984       3.557613       6.246090       8.796872    12.437498 

14.349998-11.120001    -9.296001-14.436301       -.149441 

7.080447    18.664358       5.531486      -.774811-10.019849 

_=-.9.  2 1  58  79- 15.  372703=1  0.8981  63. _6..0.8..1470._._.l.  865  176 
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23 

884 
383 

-34 

893 

-24.514 

-8.411 

-65.729 

-86 

-64 

707 

-71.165 

-37.132 

-21.906 

-43 

125 

-25 

700 

-25.960 

-79.768 

-64.414 

-20 

131 

-28, 

105 

9.316 

29.253 

22.802 

-31 

358 

-67 

287 

-65.829 

-52.663 

-43.731 

-68 

985 

-31 

.588 

-30 .670 

3.864 

49.291 

35 

633 

-3 

294 

-42.835 

-73.268 

-109.414 

-78 

131 

-24 

705 

50 -036 

21 .629 

43.703 

3  7 

562 

-23 

350 

-57.080 

-74.064 

-102.251 

-70 

001 

.-48 

801 

-32.041 

4.367 

-50.106 

-50, 

685 

-62. 

348 

-20.078 

-12.263 

16.990 

-2, 

408 

-37, 

326 

-69.461 

-38.769 

-75.415 

-95 

932 

-33 

546 

2.363 

-  .309 

25.553 

-25, 

158 

-89, 

726 

-132.581 

-87.065 

-82.452 

-87, 

562 

-32, 

649 

-2.119 

-15.695 

-25.156 

15 

275 

-43, 

180 

-16.144 

-38.715 

-42.572 

-68 

058 

-46, 

446 

-17.957 

36.834 

59.068 

12 

654 

-7, 

077 

-62.061 

-77.649 

-71.519 

-38, 

615 

2, 

908 

9.126 

-5.899 

-20.719 

-20, 

975 

-57, 

780 

-96.824 

-120.059 

-67.248 

-63, 

398 

-22, 

918 

-51.135 

-83. 308 

-97.246 

-36, 

597 

-54, 

478 

-21.782 

5.174 

55 


-15 
-62 

75 
12  7 
-85 
-46 
-22 
5 
5 
-34 
-87 
-58 
-64 

24 

2 

-74 

-39 

■17 

15 
134 
-18 

30 

-92 

-54 

-3 

•111 

16 

-76 

-4 

-9 

"-4  7 

-41 

-1 

•126 

33 


061 
015 
190 
689 
349 
064 

489 
32  6 
978 
450 
752 
047 
636 
943 
232 
886 
830 
040 
752 
747 
670 
241 
288 
336 

139" 

169 

510 

964 

928 

241 

70  5" 

574 

374 

412 

487 


-68.306 
161.561' 
137.082 
-27.003 
"-64V465 
-92.385 
1.976 

56.23  7 

■102.344 
-24.497 


-56 

120 

1 

79 

•104 

-26 

127 

-9 

19 


683' 

657 

670 

119 

9  3 1 

599 

837 

377 

131 


-2c 
-7'i 

44 

-86 
-94 
-27 

-16 
-3  0 

-e 

126 

2 

-16 

-91 

52 

106 

-45 

60 

-61 

l5y 

-2 

17 

-77 

-34 

55 

-60 

-46 

-70 

6 

-2b 

-79 

-52 

-8 

-73 

2 

6 

■18  0 

96 

1 

-52 

■123 

-68 

74 

-85 

11 

•in 

•154 

93 

'-17 

-85 

-70 

-13 

14 

-99 


248 
212 
952 
951 
480 
851 
391 
539 
582 
16  0 
001 
962 
709 
245 
986 
909 
464 
032 
198 
597 
536 
993 
430 
469 
689 
535 
592 
171 
458 
393 
764 
459 
099 
529 
590 
555 

24  9 
666 
198 
972 
5o8 
819 
.790 
675 
402 
147 
325 
536 

10  5 

184 

679 
269 
698 
695 


-52 
13 

120 
-22 
-26 
-24 
-43 

-3 
-62 
-54 
-36 

-3 
-77 
-11 
-Qi 

40 
7  7 

-92 

142 
-42 
-41 
-40 
27 

-39 
-97 

-40 
57 

-2 

160 

25 

-86 

^79 

49 

-42 

•115 

102 

-30 

-6 

■  1 3  0 

5  3 

-72 

-17 

-26 

•114 

•  14  2 

17  1 

-5 

■103 
83 

-74 
23 


999 
23  0 
438 
761 
984 
081 
113 
231 
534 
328 
201 
430 
767 
796 

OH. 
52  7 

829 
626 
958 
278 
029 
406 
344 
025 
0  49 
628 
273 
937 
966 
314 

Oil 
833 
679 
423 
872 
956 
199 
932 
642 
377 
4o7 
745 

940 
678 

517 
660 
029 
884' 

74  8 
743 

185 
042 
444 


-53 

40 

-4 

143 

-34 

13 

-42 


-66 

-51 

-4 

-69 

4 

-98 


399 
984 
951 
609 
787 
73  5 
4  91 
785 
827 
862 
761 
744 
786 
637 
-40  .009 
-5.422 
937 
301 
167 
422 
623 
325  ~ 
275 
780 
639 
297 
819 
950 
373 
148__ 
209 
666 
143 
739 
897 
564 
559  " 
746 
913 
10  2' 
725 
996 
9/5 


-9 

62 
113 
-52 
-40 
-35 
-56 
r19 
-69 

34 
-92 
-45 

49 

16 
"151 

45 
123 
-65 

51 

-96 

"-18 

90 
-38 
-45 
-42 

56 
•163 


116 

124 

138 

14 

-71 

•113 

47 

-67 

5 


-52 

.519 

100 

.987 

-86 

.361 

-103 

.687 

-11 

.830 

-74 

612 

-26 

592 

2 

972 

-41 

062 

-65 

690 

-21 

609 

-81 

795 

16 

429 

-33 

710 

-78 

607 

-24 

537 

20 

051 

27 

.441 

-117 

933 

-3 

338 

-39 

699 

-55.860 

-45.420 
-52.424 

""-80.711 
89.638 

-140 .455 


6 

52 

-88 

-40 

-7 

-14  7 

42 

-22 

-82 

1 

64 

-41 

'  "70 

-52 

67 


840 
698 
881 
967 
467 
515 
609 
882 
452 
353 
997 
33i 
481 
78q 
797 


21.848 
52.943 


-13q.18q 
-41.122 


214 
328 
623 

093 

9  9,; 

995 
348 
233 
355 


-60 

-104 

-117 

172 

-84 

-30 

-150 

60 

-48 

26 


354 
571 
663 
898 
726 
998 
796 
278 
587 
464 


56 


57 


58.799 
-125.672 

-10.142 

127.239 

-14  0.737 

-7.714 

10  2.391 

-15  6.590 

2  7.229 

109.513 

-91 .472 

-69.617 

-147  .263 

-75.Q34 
25.324 

73.410 

-41  .178 

-180 .293 

-107.435 

-47.488 

-145.102 

-208.348 
-27.990 

-   .081. 

-196.078 
-33.792 

-46.345 

-74  .610 

-90  .627 
-27.541 

101.315 
-40.834 
-83  .567 

-16. 633 
61.252 

13.933 

-11.306 

-26.598 
5.34  6 

65.555 
-13.679 

-:  .  .323 

-26.356 
61.457 
-10  .4j>_9__ 

-105.Q53 
36.022 
12.24  7 

-118.243 

-7^.382 
66  ,597 

-71.394 

-107.706 
36.478 

-94.715 
-5.365 

4  3.18  2_ . 

-17.972 
82.308 

-10 . 054 

-17.443 
-44 . 806 

16. 934 

-66.155 

-45.845 

6  7.347 

-154.324 

66.324 

-14.122 

-146.259 
-17.541 
1 2  7  .  9  0  2  __ 

-64.0  07 

-28.-32 

38.322 

-41.143 

-230 .065 

46. 633 

-80.514 

-274.252 

76.106 

-.126.611 

-151.202 

12.285 

-146.289 

-231.761 

105.223 

-182.831 

-19.209 
11.182 

25.946 

-115.527 

19  .4q2 

-1  .443 

-142.822 

104  .594 

66.845 
-93.257 

115.275 

.6  76 

-42.606 

8.420 

-103.659 
7.115 
6.336 

42.669 

-161.777 

-8.656 

10 .735 

-105.222 
-38.725 

-94.612 
-65.577 
-; 5.18  0 

-13  7.090 
4  7.538 

- 1  0  6  .  1  4  4 
8  4.066 

-142.227 
-3.573 

-20  6.472 
7.430 

-118  .915 

157.455 

-135.931 

-35.258 

-30  .132 

76.964 

-107.385 

14 .294 
93.771 

-92.70 3 
69.035 

-64.121 
-17  0  .485 

-1.164 
-16.783 
-23.966 

16. 79^ 

85.099 
-221.357 

23.428 

7.30  3 
-179.388 

-51  .658 

-33.557 

-112.51 1 _ 

64.8/4 

-189  .446 

-33.409 

-1.527 

-25.108 
1  4  4  n  5 

53,979 

-60 . 086 

67  .924 

25.383 

-17.069 
-3.261 

-52  .  294 

120.945 

-65.409 

-  1  c.  0  ,  0  0  3 

116.756 
-176.727 

-        "     -103.732 

2.172 

-  7  4  3  ft  6 

-130.025 

-45.062 
37  89  2 

-16J.820 

11.950 

-71.687 

-122.856 

27.960 

-75.149 

-21.285 

-58.232 

11.480 

91.964 
-176.205 

62.092 

47.788 

-111.564 
36.074 

3.630 
-2  0  .851 
74  ,859 
SECONDS 

-23.696 

92.519 

-90.113 

-15.757 

-5.385 
-70.490 

TIME,  14  MINUTES 

AND   37  i 
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APPENDIX  II-A 


PROGRAM  MEAN 

DIMENSION  R( 1OQ0) , E( 1 '  _ 

"HIS  PROGRAM  SIMULATES  T  H  c  B'ASYC  MODEL  WITH  MEAN  DEMAND  INPUT 
DO  1  1=1,1000 

1  R( I )=16.0 

DO  U  1=6,1000 

SUM=0.0 

'••=1-5 

DO  2  J=1,M 
TEMP1=R( I  >*(  ,8)**J 

2  SUM* SUM+ TEMPI 
TEMP1=SUM 
SUM=0.0 

DO  3  J=1,U      

TEMP~2=R  ;  I  ; 

3  SUM=SUM+TEMP2 
TEMP2=SUM 

k    E( I)=TEMP1-TEMP2 
PRINT  5 

5  FORMATt 11HMEAN  VALUES) 
PRINT  6, (E( I ), 1=6, 10001 

6  FORMATi 6X,5F10.6) 
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APPENDIX  I I -3 


PROGRAM  MEAN 

DIMENSION  H  C 100  0  )  , £< 1000  ) 
THIS  PROGRAM  SIMULATES  MlDEl  WITH  MS  AN  DEMAND  INPUT 
0  0  1  1=1.10  0  0 

1  R(I  )  =  1 

DO  4  1=6.1000  =__=_=_ 

SUM=0.0 

M=I-5 

DO  2  J  =  i,K 

TEMPl  =  RU-*-J)  -  ( 

2  SUM=SUM+T£MP1 
TEMP1=SUM 
SUM=0.0 
DO  3  J=l,5 
TEM?2  =  R(  Ijlj-JJ 
SUM«SUM*TEMK2 
T6MP2=SUM 


e>*»<  J-D 


EC  1  )*T 

PRINT 

FORMAT 

PRINT 

FORMAT 

END 

END 


EMP1 

(HHMbAN 
6,(£(I), 
(6X,5Fl3 


EMP2 


VALUfi3) 
?6,1000 
6) 


APPENDIX  II-C 


„  ^    T  PROGRAM  POISSON 

C_JQilS  PROGRAM  SIMULATES  BASIC  MODEL  WITH  RANDOM  DEMAND 

DIMENSION  T(  1000),xn  000),  IX(  1000) 

COMMON  I  RANDOM, RANDOM 

IRAND0M=3574655 
DO  'I  1=1,1000  

CALL  NORMALi  16.0, 4,0, X(  If)"7       ~~~ 

IX  < I )=X  C  I  ) 

1  X{  I)  =  IX(  I  ) 

DO  4  1  =  6,1000 

SUM"=0.0  

M=I-5 

DO  2  J=1,M 

_TEMP1=( .8)**(J-1)*X( I-U-JJ 

2  SU M= SUM  +  TEMPI 
TEMP1=SUM 
SUM=0.0 

DO  3  J=l  ,5 

TEMP2=X  (  I-J  +  IT 

3  SUM=SUM+T£MP2 
TEMP2=SUM 

h    Tt I )=TEMP1-TEMP2  

PI  NT  5  "" 

5  FORMAT! 24HRAND0M  NUMBERS  GENERATED) 
PRINT  6,  (X( I  ),  1  =  1 ,  ,000) 

6  FGRMAT{ 6X,5F10.6)  

PRINT  7 
7  FORMAT( 17HINVENT0RY  ON  HAND) 
PRINT  8,  CT( I  it  1  =  6, 1000) 

8_F0R.MA7(  6X,5F10.6) 

END 
^SUBROUTINE  RANVAR 
CONMGN  I RANDOM, RANDOM 

C0N{K1  =  67108864)       ; 

IRANDCM=IRANDCM*3125 

I F { I  RAN DOM-6 7 10  8865)2,2,1 

1  ENQ(O), LDA( I  RANDOM  ),OVI (Kl ) ,STQ{  IRANDOM) 

2  Y= I  RANDOM  

~RANDGM=Y/6710886U, D 

RETURN 

END_ 

StTTTRGUTINE    NORM  AL  (  A,  B,C  ) 

TtfMMON    IRANDCM,  RANDOM 

SUM=0.0 

DO    1     1=1,12 

CALL    RANVAR 

T~  SUM=SUM+RANDOM  " 

X=SUM-6.0 

C=X*B«-A 

RETURN 
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APPENDIX  II-D 


PROGRAM  POIaSON 
THIS  PROGRAM  SIMULATES  BASIC  MODEL  WITH  RANDOM  DEMAND 
D  I  MENS  I  ON  ~110CO),X<10GO),IX(10GO  ) 


COMMON  I  RANDOM • RANDOM 

IRA\D0M=357465b 

DO  1  1=1.10.00 


CALL  NORMAL ( 1 6 . 0 , 4  ,  0  ,  X  (  I  )  ) 
IX( I >  =  X(  I  ) 

_XJ  I  )  =  i  X  (  [  )   

DO  4  1=12,1000 
SUM  =  0  .0 

M=I-11  


DO  2  J=1,M 
TEMP1=2.2+(.8)**(J-1)*X(I-1Q-J) 

2  SUM=SUM+TEMP1      

TEM?1=SUM 
SUM  =  0-0 

DO  3  j  =  i ,  1 1  ■ '__ 

TEMP2=X( I-J+l) 

3  SUM=SUM+TEMP2 
TEMP2=SUM  ._ 

4  T(]  )s  TEMPI-  'I  EMI  2 
PRINT  5 

5  F0RMAT(24hRANDuM  NUMBERS  GENERATED) 
PRINT  6 , ( X  ( I).  1=1,1000) 

6  F0RMAT(6X,5l-  10-6) 
PRINT    7 

7  FORMAT(i7HllWENTORY    ON    HAND) 
PRINT    8,(Ki)-  ! '  =  1 2  , 1 0  0  0  ) 

8  FORMAK6X.5F10  .4)        _ 

END 

SUBROUTINE  KANvAR 

COMMON  I  RANDOM. RANDOM 

CON(K1=67108864) 

IRANDOMs (KANDOM*3l25 

IFl IRAND0M-67in8863)2,2,l 

1  ENQ(0)  *LDA*  1  RANDOM)  ,OVI(*l>, ST  KI  RANDOM) 

2  Ya I  RANDOM 
RANDOMsY/67108664  .0 
RETURN 

END 

SUBROUTINE  N0KMAL(A,8,.C) 
COMMON  IRANDOM. RANDOM 
SUM  =  0  .0 
DO  1  1=1,12 
CAL-  RANVAR 
1  SUM=SUM+KANDOM 
X=SJM-6.0 
C  =  X  *  8  *  A 
RE  1  URN 
ENJL. 

62 


\"DPTP  FlTY  TT  ' 


PROGRAM  POIbSO! 
C  THIS  PROGRAM  SIMULATES  BASIC  MODEL  WITH  RANDOM  DEMAND 
DIMENSION!  i  UOPO  ),  XCiOOO  ),  1X(  18  00  )      

COMMON  IRANDOM, RANDOM 

IRAND0M=3574655 
D  O  .  1 .  I  =  i  ,  1  0  0  0  _____ 

CALL  NORMAL ( 100  .  0 , 10  .  0 # X ( I > ) 

I  X  (  I  )  =  X  (  I  ) 

i_x(i)  =  ixd) 

DO  4  1=6,1000 
SUM  =  0  -0 

M=I-5 

DO  2  J=_,M 

J  )** ( J-l  )*X(  I -4 
EM  Pi 


2  J  =  l 
T£MP1=( . 
SUM=SUM+ 
TEMP1=SUM 
SUM=0-0 
DO  3  Jr_,5 


J) 


6 


TEMP2=X( I-J+l) 

SUM=SUM+TEMP2 

TEMP2=SUM 

T (  I  )  =  TEMP1- i EMP2 

PRINT  5 

F0RMAT(24HRAND0M  NUMBERS  GENERA  TED L 

PRINT  6,  (X( I  ) ,  1=1,1000) 

FORMai <6X,5F10 .3) 

PRINT  7 


7  FORMAT (17H INVENTORY  0\  HAND) 


8,(1(1), 
.<6X,5»F10 


=6.1000 ) 
4) 


AR 
RANDOM 


PRINT 
8_  FORMA 
END 

SUBROUTINE  KAN 
COMMON  IRANDOM 
CON(K1=67108864) 
IRAND0M=IRAND0«v*3l25 
LFJ  I  RANDOM- 671  0  8863)  2,  2,1 ■ 

1  ENQ(0 ) ,  LDA( IRANDOM) , DVl(Kl) , STQ( IRANDOM) 

2  Ys I  RANDOM 

RAND0M=Y/671q8864.q   

RETURN 

END 
SUBROUTINE  NORj^ALU, 

COMMON  IRAN uOM, RANDOM 

SUM=0 . 0 

DO  1  1=1,12      

CALL  RANVAR 
1  SUM=SUM+RANDOM 

X=SUM-6.0  

C=X*B-A 

RETURN 

END  
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,:o:  I--1* 


PROGRAM  POIbSOM  nnvi  nCMAMn 

THIS  PROGRAM  SIMULATES  BASIC  MODEL  WITH  RANDOM  DEMAND 
DJ.ME.NJS_I.0R  1  HOOO),X(1000  L_  IX  (100  0) 
COMMON  I  RANDOM.  RANDOM 
IRAND0M=3574655 

DO  1  1=1,100  0 

Z=l 

W  =  16.  0+Z 
U=SQRiF(w) 

CALL  NORMAL (W.U,X< I > > 
I  X  < I >  =x  C 1  ) 
1  x  (  I  )  =  I  X  (  I  ) 


DO  4  1=6,1000 

SUM =0-0 

M  =  I  -  o 


DO  2  J=1,M 
TEMPl=(.8)**U-l>*X<I-*-J> 

2_SUM.=  SUM  +  TEMPI 
TEMP1=SUM 
SUM  =  0  .0 
no  3  J  =  l»5 
TEMP2sX( I-J+l) 

3  SUM=SUM+TEMP2 

TPMP2sSUM 

4  T  (I)  =  iEMPl-  IEMF'2 
PRINT  5 

5  FOR ■:•,,-..  ,  (  2  4  h  P  A  h:  D ' ?.__  NUMBERS  GENERATED) 

',INT  6,  fx<  1  )  ,  1=1,1000) 

6  FORMAT <6X,5F 10  3) 
PRINT  7 

7  F0RMAK17HINVENT0RY  ON  HAND) 
PRINT  8, U ( 1 ) , 1=6,1000) 

8  PORMAT(6X,i?K10.3) 

END 

SUBROUTINE  KAN'-'AR 

COMMON  I  RANDOM .RANDOM 

CON(K1=67108864) 

IRAND0M=lRANDOM*3l25 

IF( IRANDOM- 67108863)2, <,1 

1  ENQ(0),LDA(IPA,DOM),DVKKl),ST.CiRANDOM) 

2  YsIRANDOM 

RANDOM  =  Y/67108864,_Q 

RETURN 

END 

SUBROUTINE  NORMAL (A, 3, C) 
^COMMON  IRANDUM. RANDOM 
SUM=0 .0 
DO  1  1=1,12 
CALL  RANVAK 
SUM=SUM*RANDOM 
X.sS.UM-6.  0 

C=X*B*A 

RETURN 

END 
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APPENDIX  III-A    - 
RANDOM    NUMBERS    GENERATED 

12-'-''05?92  4S°250000  15.000000  17.000000  8.000000 

92-2S2R8g  Jg-oopooo  15.000000  12.000000  19.000000 

20-OoOOuO  18.000000  18.000000  20.000000  14.000000 

]9*?999S9  2o.g            13.000000  14.000000  12.000000 

11.000000  22.000000  12.000000  13.000000  16.000000 

J9-999k99  t§ ."oooooo  13.000000  23.000000  22.000000 

U-222222  8.990000  11. 000000  12.000000  is.oooooo- 

J7.000D00  15.000000    18.000000    21.000000    18.000000 

J?.«000300  16.000000       9.000000    13.000000    20.000000 

16.000000  11.000000    15.000000       7.000000    15.000000 

2d. 000000  18.000000    19.000000    13.000000    19.000000 

13.000000  12.000000    15.000000    22.000000    18.000000 

20.000000  11.000000    24.000000    13.000000     14.000000 

18.000000  12.000000    14.000000    13.000000    15.000000 

2^.0000  00  11.00  0000    16,0000  00    11.00  0000    19.000000 

12.0000  00  18.000000    22.000000    11.000000     16.000000 

19.000000  12.000000     13.000000    18.000000    11.000000 

T2. 000000  12.000000    10.000000    21.000000    14.000000 

11.000000  19.000000    20.000000    23.000000    13.000000 

15.0000  00  13.000000       6.000000    16.000000    24.000000 

23.000000  15.0                  12. 0 OC     DO    24.000000    13.000000 

18.000000  19.000000    20.000000    12.000000    16.000000 
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